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Abstract: Saltation is the major particle movement type in wind erosion process. Saltating sand grains can 

rebound up to tens of times larger in length and height over hard surface (such as gravel surface) than over 

loose sand surface. Gravels usually have different faces, causing distinct response of the impacting grains, 

but the effects of the grain and gravel-surface contact angle on grain rebound are not yet well quantified. 

We performed full-range controlled experiments of grain saltation using different contact angles, grain sizes 

and impact speeds in still air, to show that contact angle increases the height of representative saltation path 

but decreases particle travel length. The results were compared with outputs from the COMprehensive 

numerical model of SALTation (COMSALT). Large saltation height of 4.8 m and length of 9.0 m were 

recorded. The maximum and representative saltation height over the gravel surface were found to be about 

4.9 times and 12.8 times those over the loose sandy surface, respectively. The maximum saltation length 

may be reduced by 58% and the representative saltation height may be increased by 77% as contact angle 

increases from 20° to 40°. We further showed that the collision inertia contributes 60% of the saltation 

length, and wind contributes to the other 40%. These quantitative findings have important implications for 

modeling saltation trajectory over gravel surface. 
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1 Introduction 


Saltation is the main mode of aeolian sand transport during wind erosion process in planetary 
environments (Shao et al., 1993; Ayoub et al., 2014; Burr et al., 2015). Saltation trajectory may be 
one of the main factors influencing the initial size of aeolian dune (Pahtz et al., 2015; Lammel and 
Kroy, 2017) and deserves the most research consideration (Li et al., 2020). The characteristic 
length and height of saltation trajectories over loose sand surfaces are at magnitude order of 10° to 
10° m (Williams and Leer, 1995; Hagen et al., 2010), and the rebound of particles over massive 
particles in a heterogeneous surface travels higher than that in a homogeneous surface with the 
same grain size (Bagnold, 1973). In the blown sand condition, smaller particles saltate higher or 
suspend in air, because the drag from wind varies as the square of particle diameter but the weight 
varies as the cube of it (Lorenz and Zimbelman, 2014). 
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However, over immobile hard surfaces, such as gravel (stony) Gobi desert that covers an area of 
1.3x10° km? in Central Asia, saltation can be remarkably different, with much longer and higher 
trajectories, indicating the effect of particle rebound with hard surface. This fact was known long 
ago, as Bagnold (1941) noted that "the sand movement over a rough immobile or partially immobile 
surface differs in this important respect from that over a surface of dune sand...", because in this 
case grains impact the surface in bouncing rather than splashing and more saltation kinetic energy 
is retained. Over a granular sandy surface, most of the grains are ejected with low energy (Ammi 
et al., 2009), but the majority particles can be rebounded by the rigid gravel surface. As a result, 
sand flux over the gravel surface under extremely strong wind conditions (gust wind speed 
exceeding 40 m/s) may be the largest in the layer tens of centimeters above the ground surface (Fig. 
1), and sand grains jumping up to 10 m high have been observed (Yin, 1989; Cheng et al., 2016; 
Tan et al., 2020). 
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Fig. 1 Evidence of large scale saltation over gravel surface. (a), serious sand abrasion of a wood pole on windward 
side from about 0.1 m to more than 1.5 m in a Gobi area of Xinjiang, but the bottom part below 0.1 m is less damaged; 
(b), abrasion of a concrete bridge pier (indicated by the ellipse) shows the strongest abrasion at about 0.8 m high; (c), 
sand deposited behind a 3.3 m high wall as wind breaker at a section of the Lanzhou-Xinjiang high-speed railway; 
and (d), grain-size distribution of five samples from the site shown in (c). At the three sites of (a)-(c), wind is almost 
unidirectional, with sand blowing toward the abrasion side in (a) and (b) but another side of the wall in (c). 


Most studies on sand movement over gravel surfaces use the same approach as that over sandy 
surfaces, e.g., examining aerodynamic roughness length and wind velocity profile (Dong et al., 
2002a, b; Zhang et al., 2008). A turning point was observed in sand concentration profile due to 
sand-gravel collision, which differs from the concentration profiles over sand surfaces (Zhang et 
al., 2008). Particles with different sizes may also respond diversely during the collision with 
surface. Although Bagnold (1941) noted that the size and spacing of the pebbles affect particle 
movement, few studies focused on how pebble characteristics influence grain movement. Hardly 
any such studies were carried out at large scale conditions due to practical difficulties; for example, 
saltation height could exceed the wind tunnel upper boundary, so that the actual height cannot be 
measured. In field conditions, grain size and gravel characteristics are still hard to be measured and 
parameterized, because it is not feasible to separate the effects of gravels with different shapes, 
sizes and surface angles on saltation particles. As a result, new experimental method needs to be 
developed to study the effect of gravel surface on particle saltation. 

Suppose a group of grains with the same properties and same trajectory (i.e., with the same 
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impact velocity and angle) impact a gravel surface. Then, the rebound would substantially affect 
the following trajectory of the grains. A natural gravel surface can be considered to be composed 
of many surfaces with different reflection angles. If the effect of the contact angle between the 
impacting grain and the gravel surface is known, the complicated trajectories of sand grains 
observed in wind-tunnel and/or field experiments can be decoded by considering the probability 
distribution of the reflection angles of the gravel surface. We can then better predict grain 
trajectories and sand transport flux. This can increase our knowledge of wavelength of incipient 
ripples or dunes in areas of hard surfaces, which is of great importance for our understanding of 
gravel surface formation and evolution. 

To obtain the trajectory of saltation sand that may not be able to get through traditional wind 
tunnel experiment and field measurement, we designed a set of equipment to obtain the full-scale 
saltation range of sand particles, and performed controlled experiments to measure the 
concentration of saltation particles and the length of grain saltation at sequential locations. We 
obtained details of saltation trajectories in still air to show how contact angle of gravel surface 
influences sand saltation trajectories for different grain sizes and impact speeds. By analyzing the 
characteristics of the trajectories, we distinguished the contributions of grain inertia and fluid 
momentum to sand saltation. The results were then compared with simulation for sandy surface 
from the COMprehensive numerical model of SALTation (COMSALT), a physically based 
saltation model. which is able to switch off the effect of turbulent air and generate comparable 
results against our tests (Kok and Renno, 2009). 


2 Materials and methods 


2.1 Blasting equipment 


We designed an equipment that can control the blasting speed of air and contacting angle, and 
collect the impacted particles at various horizontal and vertical locations (Fig. 2a). An air-compressor 
was used as the power section to create high-speed air flow; a pressure valve, to control air pressure; 
an electromagnetic valve with a timer, to control blasting time; a sand-blasting gun with a round 
opening of 0.5 cm in diameter, to feed and eject sand particles; a hard smooth slate with adjustable 
angles on a platform, to simulate the gravel surface; a series of square plastic bins (10.0 cm in length 
and 30.0 cm in depth) along the blasting center, to collect the horizontally settled particles (Fig. 2b); 
and a stack of removable steel bins (each 5.5 cm in height, 5.0 cm in width, 20.0 cm in length, 25° 
uptilt, with a polyporous air outlet at the end) on a portable frame, to collect the saltation sand profiles 
at different distances (Fig. 2c). 

The air speed was calibrated at 15.0 cm (the location of proposed contacting center) from the 
opening of the sand-blasting gun. A pilot tube was used to measure the wind speed change with 
controlled air pressure and build their relationship. The pressure was increased from 1.0x10° Pa at 
0.5x10° Pa interval with three replications. A good linear relationship was determined (Fig. 2d); but 
above 3.5x10° Pa pressure with about 56.0 m/s air speed, the relationship tended to degrade. In fact, 
the pilot tube failed to reflect any air speed change with higher than 4.0 10° Pa pressure, because the 
measured static pressure also increased in the violent airflow. So, only four reference air speeds from 
20.0 to 50.0 m/s with interval of 10.0 m/s were used. The selected speeds may be apparently higher 
than those in blown sand experiment, especially in wind tunnel, but can be frequently visited at Gobi 
areas. The particles travelled a short distance (15.0 cm) in the same direction together with the high- 
pressure jet air stream, so we can assume they followed the speed of the air before collision, with 
little damping during the transient travelling period (<0.0075 s at 20.0 m/s air speed). 

The vertical bins were arranged with their openings on a vertical plane, and the lower edge of the 
lowest bin was set at the same height with the contact center. The dimension of the contacting slate 
was decided to make sure that no particles can directly fly over but have to collide with the slate at 
the three angle conditions (Section 2.3). Trial tests were carried out before each speed, angle and 
particle size experiment to decide the number of needed bins at the two directions. For each scenario 
with a certain contact angle, speed and particle size, the horizontal test was performed first, and then 
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the vertical one with the frame stood from the margin of the slate (0.0 cm) and every 30.0 cm moving 
back to obtain the trajectories. The bins were taken out starting from the top ones, and sand particles 
in each bin were weighted by a ten thousandth balance. If particles were found less than the maximum 
resolution of the balance, their numbers were counted and their weights were estimated according to 
the density of quartz and average volume (assume sphere shape) of the size group. Three replications 
were performed for each test scenario. 

All the tests and weighting works were carried out in indoor environment with still air. During 
each test, a sheet of cloth was laid on the floor to recycle the rest of blasted sands and also to prevent 
secondary saltation. 


2.2 Contact angle at gravel surface 


In the Gobi area, gravels on the underlying surface have different coverage, sizes, orientation, and 
roundness. The shapes of the gravels can be very irregular, and each gravel may have several faces. 
The inclination angle (0 for horizontal) of each face that may encounter a saltation sand defines the 
contact angle 6; then, the rebound liftoff angle (a) will be the sum of 26 and £, where £ is the 
blasting angle of particle (downward is positive). In the used equipment, J is set to 0 (horizontal 
stream of sand) for simplification, so a=20. 

We randomly picked over 300 gravels from the Gobi in eastern Kumtag Desert, Northwest China 
and manually measured their faces. Beside the underside face that touches the ground, each gravel 
usually has 3—4 faces that may encounter saltation sands. The face inclination angles (0° is at level 
with ground, and 90° is vertical) distribute over a large range from almost 0° to even over 90° (Fig. 
2e). But for all the steep gravel face angles that are larger than 50°, the grains would rebound 
vertically or even backwardly upwards so that their forward velocity and next travelling paths 
cannot be measured. For the gentle angles smaller than 10°, the saltation condition would be similar 
to that over a level surface, and the grains may not be able to develop a full trajectory before hitting 
the ground. As a result, we only considered gravel inclination angle of 20°, 30° and 40°, which 
were also the most distributed measured angles and resulted in the a of 40°, 60° and 80°, 
respectively, assuming elastic collision. The air out of the nozzle may diffuse at a degree about 13° 
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Fig. 2. Equipments used in the experiments. (a), a sketch of the test equipment with the power section, control 
panel, contact slate with its contact angle (9), a horizontal nozzle to simulate an impact angle (£), the liftoff angle (a) 
of rebound particles, and the horizontal and vertical sand-catch bins (the inset shows the relationships of 0, a and 2 in 
details); (b), horizontal bins to collect sands on the floor; and (c), removing of vertical bins; (d), tested relationship 
between used pressure and air speed, in which v is impact speed and p is air pressure; (e), distribution of counted 
gravel inclination angles. 
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and cause errors (Reichardt, 1941), but sand particles have much larger inertia and result in less 
dispersed a than the air. 


2.3 Sand 


Aeolian sands were taken from the Mingsha Megadune to the east of the Kumtag Desert, Northwest 
China. These quartz sands were screened to form three size groups: 0.10—0.25, 0.25—0.50 and 0.50- 
1.00 mm as fine, medium and coarse sand, respectively. 

In total, 36 groups of experiments were carried out. Each involved removal and reload of tens of 
horizontal bins for one time and vertical bins for 6 (the shortest measured distance, 1.8 m) to 27 
(the longest measured distance, 8.1 m) times. 


2.4 COMSALT simulation 


The experimental results on saltation with different 9 were compared with simulations from 
COMSALT (Kok and Renno, 2009). This model simulates particle saltation due to gravity, fluid 
drag, fluid shear, turbulent air, and a wide range of other processes. It does not consider the collision 
of saltation particles with each other; and it can turn off the effect of turbulent air, which makes it 
similar to our experiments in still air. Many fluid forces like the Saffman force from the shearing 
flow and the Magnus force from particle rotation were experienced in both the experiments and 
COMSALT. The model reports two-dimensional steady-state results, meaning that interactions of 
particles with surface and wind are averaged over many iterations so that a dynamic balance is 
reached. 

COMSALT assumes a flat surface and estimates a mean rebound angle a of 40°, which fits our 
6=20° experiments. As a result, the comparison will show the difference of saltation from gravel 
surface (our test) against flat surface (COMSALT simulation). However, it is worth noting that 
other ejection angle ranges were also recommended, for example 60° (Ammi et al., 2009). 

COMSALT uses the log-law of wind profile to estimate the shear velocity of clear wind using 


u z 
u=—In—, 1 
"i (1) 
where u is air speed (m/s); u” is shear velocity; z is height (m); zo¥D/30 is surface roughness based 
on particle size D (m); and k is van-Karman constant and set to 0.4. Details about particle trajectory, 


collision process and readjustment of wind profile can be found in the document of this model (Kok 
and Renno, 2009). 


3 Results 


3.1 Distribution of saltation sands 


Since sand flux was continuously measured, we can plot the concentration contours of the saltation 
sands, which also represent the probability density of sand grains for given size, impact speed and 
6 conditions. For all 6 conditions, sand grains tend to concentrate in a narrow range of trajectories, 
indicating that there exists a dominant saltation trajectory where particle concentration can be 
hundreds of times larger than those at the margin area, but the scatter range of particle saltation 
(light gray areas in Figs. 3—5) is obviously large. This is because saltation process is highly 
stochastic, and the particles are unlikely to transport in one path even with the same internal and 
environmental factors. 

As 0 increases, particles move more upward but settle at shorter distances, resulting in larger 
saltation range after collision and more upright (higher) but narrow (shorter) profiles, as in Figure 
3a at 20° vs Figure 3c at 40°, for example. Also, the peak concentration of particle cloud is reduced 
more quickly with 6 although the total flux remains the same with the same air speed (same inject 
speed and time duration), e.g., from 88.2 g/min in Figure 3g at 20° to 7.1 g/min in Figure 3i at 40°. 
This is because a shallow horizontal 0 gives small a, so that particle streams tend to concentrate in 
a narrow path with small 8. 

As the impact speed increases, both saltation height and length increase. The maximum saltation 
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height (Hm) and maximum saltation length (Lm) at 50 m/s impact speed for the three size groups are 
2.8 and 2.7 times that at 20 m/s wind. Larger particle sizes rebound higher and saltate longer, so 
that the saltation space follows the sequence of coarse sand (Fig. 5)>medium sand (Fig. 4)>fine 
sand (Fig. 3). We find that coarse sands with 0.50-1.00 mm in size rebound the farthest when 
6=20°, while the highest is at 06=40° under impact speed of 50 m/s. 
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Fig.3 Concentration contour maps of the saltation space for fine sand. Each column shows different impact speed 
(v) with the same contact angle 0, and each row shows the effect of 0 under the same impact speed. Large light gray 
areas in each penal indicate low possibility to encounter particles, and most of the particles concentrate at a small 
range (the representative trajectory). Log scale is used to show the changes of region with low concentrations. Some 
isolated light patches exist, because zero concentration was measured below or before the patches. 
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Fig. 4 Concentration contour maps of the saltation space for medium sand 
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Fig. 5 Concentration contour maps of the saltation space for coarse sand 


3.2 Representative trajectory and maximum height profile 


The height of maximum saltation sand concentration (He) at each horizontal location was extracted 
to form the representative trajectory (characteristic grain path) of sand particle with different size, 
contact angle and impact speed (Fig. 6). The outer boundaries of Figures 3—5 were extracted to create 
a group of Hm profiles that represent the upper limit positions of grains along the travelling paths. 
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Fig. 6 Profiles of maximum concentration (Hc) and maximum position (Hm). Each column shows different 6, and 
each row shows different particle sizes. Hm is apparently higher than Hc, meaning that particles may rise to a much 
higher location. 
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The Hn profiles bounded by the highest particles are remarkably larger than the He profiles that 
represent the average particle features. Peak heights on the Hm profiles are about 3.5, 2.7 and 2.3 
times of those on the H. profiles when 6 equals to 20°, 30° and 40°, respectively. So, if we only 
consider a representative height on the theoretical trajectory, we may largely underestimate the 
potentials of sand saltation height. The results also show that representative saltation paths have 
elliptic shape and long tails with the downward settling angle smaller than the liftoff angle 
(Bagnold, 1941; Ho et al., 2014), and the increase of H. height with wind speed (Zheng et al., 
2004). 


3.3 Characteristic values of the trajectories 


The peak Lm, He and Hm values in the saltation spaces were extracted to show their relationships 
with 0 (Fig. 7). On average, Lm increases with particle size and decreases with 0 for a given flow 
speed (Fig. 7a—c), reduced by 12% as 0 increases from 20° to 30° and by another 24% as 0 increases 
from 30° to 40°. He increases with particle size as well as with 0, except for the case of impact 
speed 50 m/s with 6=30° (Fig. 7d—f). Hm also increases with particle size and speed; however, its 
relationship with 8 is not strong (Fig. 7g—1). 

Within a size and impact speed group, 0 leads to 58% (40%—98%) variance of Lm, 77% (29%— 
97%) variance of He, but only 21% (8%—39%) change of Hm. This indicates that Hm is less impact 
by contact angle, so that particles in all size groups can reach high places after collision. 
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Fig. 7 Characteristic trajectory values for three grain sizes (one column for each size) under different impact 
speeds and contact angles. Lm, maximum saltation length. 


The changes of Lm, He and Hm values with impact speed are generally monotonic for any particle 
size and a, so that linear relationships can be established with v. The regression coefficients, a and 
b, together with their R? are given in Table 1 to show these changes quantitatively. For instance, the 
reduced values of coefficient a with increasing a for Lm indicate that trajectory length is reduced at 
a higher contact surface angle. 
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Table 1 Linear regression coefficients for Lm, He and Hm with impact speed v at different particle size and contact 
angle 


y=axvtb 
Particle size 0 La H. Hn 
a b R? P a b R? P a b R? P 
20° 13.40 -74.00 0.92 0.04 1.87 -13.20 0.93 0.03 4.84 1.10 0.85 0.07 


Fine 30° 9.90 -24.00 0.77 0.12 3.14 -34.10 0.83 0.09 6.05 -24.75 0.69 0.17 
40° 8.10 -36.00 0.88 0.06 1.38 24.75 0.95 0.02 6.88 —49.50 0.87 0.07 

20° 17.90 -124.00 0.98 0.01 2.26 -17.05 0.99 0.00 10.34 -100.65 0.98 0.01 

Medium 30° 11.80 -3.00 0.84 0.08 2.37 -165 0.98 0.01 6.99 7.15 0.91 0.04 
40° 13.40 -134.00 0.97 0.01 4.24 -34.10 1.00 0.00 9.19 -50.60 0.98 0.01 

20° 16.50 85.00 0.99 0.00 2.08 15.50 0.94 0.03 8.75 -32.45 0.97 0.01 

Coarse 30° 14.40 111.00 0.96 0.02 2.20 57.75 0.56 0.25 9.19 9.90 1.00 0.00 
40° 12.90 -14.00 0.99 0.00 3.57 35.75 0.83 0.09 10.01 -28.60 1.00 0.00 


Note: 0, contact angle; a, regression coefficient; v, impact speed; b, regression constant; L,,, maximum saltation length; H., maximum 
saltation sand concentration; H„, maximum saltation height. 


3.4 Relationships of a with L,,/H, and Lm/Hm 


A good linear relationship between L,,/H- and a was found in a previous test of blown sand in wind 
tunnel experiments, which showed that the representative saltation length can be 10.0 times of the 
height (Ling and Wu, 1980), at least when 30°<a<90° (Zou et al., 1999), which covers the current 
a range. The average Lm/He value of four speeds at three a conditions in the experiments is 40%+5% 
less than their blown sand counterpart, but they also follow a linear pattern (Fig. 8). This indicates 
Lm under a blown-sand condition in the field may be 1.7 times (1.0/0.6) larger than the current 
measured value. 

Because Hn is clearly larger than He, the ratio of Lm/Hm is smaller compared to the ratio of Ln/He, 
even after adjusting Lm by 1.7 times. The scatters of Lm/He and Lm/Hm generally distribute in 
separate areas (Fig. 8). However, due to the fact that Hm was hard to obtain in field measurement 
or wind tunnel experiment, this Lm/Hm vs a relationship is less valuable compared to the one of 
Ln! He vs a, but Lm/Hm also follows a linear changing pattern with a. 

This relationships of a with Lm/He and Lm/Hm are generally independent from wind as shown in 
Table 1 and Figure 8, which indicates that at least in the tested liftoff angle range, the shapes of the 
grain trajectories under different winds remain similar. 
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Fig. 8 Relationship between saltation length/height (Lm/Hc and Lm/Hm) and liftoff angle (a) of the test results in 


still air. Unfilled symbols are those adjusted by 1.7 times to match blown sand results in wind tunnel experiments 
by Ling and Wu (1980). 
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3.5 Comparison with COMSALT-simulated saltation over flat sandy surface 


The simulation results of COMSALT agree well with the 6=20° tests on Lm for fine sand at lower 
impact speeds (<40 m/s), but the agreement is reduced for larger particle sizes (discrepancy about 
60% for medium sand and 100% for coarse sand) although the same changing trends exist (Fig. 
9a). A possible explanation for the discrepancy is that COMSALT may overestimate L,, for large 
particles with large shear velocity (u is from 0.63 to 1.78 m/s), but may also due to the fact that 
our tests failed to collect all the saltation particles at the downstream ends where mass flux was 
very low. 

However, COMSALT clearly underestimates the saltation height over hard surface for both the 
maximum one (Hm) (Fig. 9b) and the representative one (He) (Fig. 9c). Here, we use the height of 
Zso (the saltation layer height below which half of the horizontal mass flux occurs, and the typical 
height of majority saltation particles) from COMSALT output as counterpart of our He. The Hn test 
results are 4.9 times of the COMSALT simulation results, and the He results are 12.8 times. The 
differences are similar to those between the experiments referenced by COMSALT over flat sandy 
surface with centimeter-level saltation height and the hyper saltation over Gobi areas with meter- 
level height as shown in Figure 1. This indicates that the saltation height over the gravel surface is 
clearly higher than that over the sandy surface despite similar rebound angle, and additional 
parameters are needed to reflect the saltation over the gravel surface. 


20.0 
(a) —s—COMSALT 


15.0 —o— Test 


10.0 


L,, (m) 


0.0 as 
4.0 (b) 
3.0 

& 20 

x 10 Ce 


H, (m) 


= 20 30 40 50 20 30 40 50 20 30 40 50 
Fine Medium Coarse 
Impact speed (m/s) 


Fig. 9 Comparison between the test and COMSALT-simulated Lm (a), Hm (b) and He (c) in still air for three 
particle size groups. The test is found to be conservative in reflecting saltation length (a), but clearly shows that 
current saltation model for flat surface conditions underestimate the saltation height over hard contact surface like 
the Gobi. The differences can be 4.9 times for Hm and 12.8 times for He. 


4 Discussion 


The previous section clearly showed that the effect of grain size and impact speed on saltation, and 
that the height and length of trajectory profile enlarge with the size and speed variables. In this 
section we focus on the effect of contact angle and the contribution of particle inertia on saltation, 
which are not usually available through common saltation research method. 
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4.1 Effect of contact angle on saltation 


Direct measurement of full-scale grain trajectory over gravel surface has not been performed 
before. Sand transport above a certain height, i.e., at 0.6 m, is usually suggested to be neglected 
due to its trace rate (Dong and Qian, 2007). The average saltation height is often expressed by grain 
size and wind velocity. It is true that saltation density reduces quickly with height at the space 
above the representative trajectory, as shown in Figures 3—5. Both Lm and He are affected by the 
contact angle and rebound angle. For example, here Lm decreases by about 58% when contact angle 
changes from 20° to 40°, while H. increases by about 77% with similar contact angle changes. 
They all increase with impact speed in the quiescent air tests as those in windy conditions (Ho et 
al., 2011). 

The decrease of L with @ in the still air condition can be easily imaged, and the results prove the 
reliability of the used equipment and method. Then, the results reveal that Hm is less sensitive to 
contact angle, so almost all the particle size groups can reach a high space under the same wind 
condition (more than 2 or 3 times of He), although the density or possibility is hundreds of times 
lower. The fine sand can saltate over 3.0 m high, and the coarse sand, almost 5.0 m high in the test 
scenarios, which are more than 10.0 times of the usually obtained saltation heights over flat 
surfaces. Given that wind may be more violent than the test condition and the area is large with 
unlimited sediment source, it is not surprising that aeolian sand can deposit on the trails after a 3.3- 
m high windbreak wall in the Gobi area in Northwest China, and that all the size groups were found 
as other normal aeolian sands (Fig. 1d). This explains why the sand control measure that usually 
works for sandy desert fails in some Gobi areas with strong wind. 

These findings are different from the common trajectories of blown sands, where those rebound 
with larger angles may rise higher and travel longer due to the more transferred momentum. So, 
the contribution to the whole particle trajectory can be divided into two parts, that from the inertia 
after collision and that from the wind energy. Now, we are able to separate these two contributions 
quantitatively. 


4.2 Contribution of particle inertia to trajectory 


The blown sands are supported by both the inertia after collision with surface and the energy from 
wind. To date, the available models on saltation usually adopted forces at the horizontal and vertical 
directions to solve the trajectory, and the effects of inertia and wind energy were mixed together. 
In the current test, all the momentum of sands was from the inertia, and the horizontal velocity 
increase by wind was not included in the test under the still air condition, so we can identify the 
influences of inertia and wind energy in shaping the saltation trajectory. 

The linear relationships of L,,/H. and a in Figure 8 are independent from impact speed, which 
reflect the saltation mechanism. The gap of travel distance between blown sand and still-air sand 
is due to the momentum transfer from the wind. There was no wind force to create any acceleration 
in the still-air case, but particles still travelled long distances due to the inertia and followed the 
linear relationship between trajectory characteristic value and a. Based on the offset of Lm/He values 
in Figure 8, we can deduce that during the hyper saltation process after collision with a hard surface, 
60% of the travelling path is accomplished by the inertia after the collision, and the wind energy 
contributes to the other 40%. The total Lm for the tested sands in nature wind would be about 1.7 
times of the measured size. 

The find that inertia contribute to 60% of the saltation trajectory can be used to validate the 
available trajectory models in the literature, like the distribution of lift grains and the contribution 
of the inertia. The models can consider the linear relationship of Lm/ He and a by first reflecting the 
inertia after collision, and then adding the effect of wind. In this way, the models can simulate the 
saltation process in quiescent air (when wind cases) or strong wind over hard surfaces, as well as 
the dynamic change of saltation trajectory when wind speed changes. 

The additional increase of saltation height and length of sand grains by Magnus force from the 
spin of a particle could range from 10% to about 50% (Zou et al., 2007). As a result, the test result 
may still be conservative in projecting the saltation length of sand grains over the gravel surface. 
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5 Conclusions 


Full-scale indoor experiments are performed on the saltation of sand particles over hard surface 
with different contact angles in still air. We showed that all the size groups of sands can rebound to 
several meters high, which lead to the surprisingly strong abrasion at high positions and deposition 
after high sand barriers in the Gobi area, and larger particles tend to rebound higher over longer 
distance. When the contact angles (0) are 40° and 20°, the measured saltation height and length 
reached to 4.8 and 9.0 m, respectively, with coarse sand at 50 m/s test impact speed (Fig. 7). The 
maximum saltation length (Lm) may reduce by 58% and the representative saltation height (He) 
increase by 77% when contact angle of surface changes from 20° to 40°. The collision inertia with 
hard surface contributes 60% of the saltation length, and wind contributes 40%, so the actual 
saltation length and height in the field condition are about 1.7 times larger than the measured values. 
Linear relationships are built for estimating the trajectory characteristics of saltating sands with 
different impact speed and contact angles over gravels in still air (Table 1). 

The results were compared with simulations from COMSALT, which showed that the saltation 
over hard surface is clearly different from that over a sandy surface, given all other similar factors. 
The maximum saltation height (Hm) and H: over the hard surface might be 4.9 times and 12.8 times 
of those over the sandy surface. These findings are not feasible to get in current wind tunnel 
experiments or field measurements, and can be used as the first step to understand the effect of 
surface element (gravel shape, orientation and shape, etc.) on sand saltation. New parameters and 
modules can be added to current saltation models based on the findings in this paper, and a 
physically based parameterization of particle movement over different surfaces (sandy and gravel) 
can be developed. Also, some advanced testing technology, such as particle tracking velocimetry 
(Jiang et al., 2017), might be able to work together with the sand-blasting method and measure the 
velocity field of saltation sand. 


Acknowledgements 


This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences 
(XD23060201), the Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering 
Guangdong Laboratory (Guangzhou) (GML2019ZD0601), the National Natural Science Foundation of China 
(42071014), and the Excellent Youth Innovation Promotion Association of the Chinese Academy of Sciences 
(Y202085). We thank Professor Ralph LORENZ from Johns Hopkins University's Applied Physics Laboratory and 
Professor SHAO Yaping from University of Cologne for useful discussions and comments. 


References 


Ammi M, Oger L, Beladjine D, et al. 2009. Three-dimensional analysis of the collision process of a bead on a granular packing. 
Physical Review, E79: 021305, doi: 10.1103/PhysRevE.79.021305. 

Ayoub F, Avouac J P, Newman C E, et al. 2014. Threshold for sand mobility on Mars calibrated from seasonal variations of sand 
flux. Nature Communications, 5(1): 5096, doi: 10.1038/ncomms6096. 

Bagnold R A. 1941. The Physics of Blown Sand and Desert Dunes. London: Chapman & Hall, 72. 

Bagnold R A. 1973. The nature of saltation and of 'bed-load' transport in water. Proceedings of the Royal Society of London. 
Mathematical and Physical Sciences, 332: 473-504. 

Burr D M, Bridges N T, Marshall J R, et al. 2015. Higher-than-predicted saltation threshold wind speeds on Titan. Nature, 517: 
60-63. 

Cheng J, Lei J, Li S, et al. 2016. Disturbance of the inclined inserting-type sand fence to wind—sand flow fields and its sand 
control characteristics. Aeolian Research, 21: 139-150. 

Dong Z, Liu X P, Wang X M. 2002a. Aerodynamic roughness of gravel surfaces. Geomorphology, 43(1—2): 17-31. 

Dong Z, Wang H T, Liu X P, et al. 2002b. Velocity profile of a sand cloud blowing over a gravel surface. Geomorphology, 45(3— 
4): 277-289. 

Dong Z, Qian G Q. 2007. Characterizing the height profile of the flux of wind-eroded sediment. Environmental Geology, 51: 
835-845. 


Hagen L J, van Pelt S, Sharratt B. 2010. Estimating the saltation and suspension components from field wind erosion. Aeolian 


LIU Benli et al.: Large scale sand saltation over hard surface: a controlled experiment in still air 


Research, 1(3—4): 147-153. 

Ho T D, Valance A, Dupont P, et al. 2011. Scaling laws in aeolian sand transport. Physical Review Letters, 106(9): 094501, doi: 
10.1103/PhysRevLett.106.094501. 

Ho T D, Valance A, Dupont P, et al. 2014. Aeolian sand transport: Length and height distributions of saltation trajectories. Aeolian 
Research, 12: 65-74. 

Jiang C, Dong Z, Wang X. 2017. An improved particle tracking velocimetry (PTV) technique to evaluate the velocity field of 
saltating particles. Journal of Arid Land, 9(5): 727-742. 

Kok J F, Renno N O. 2009. A comprehensive numerical model of steady state saltation (COMSALT). Journal of Geophysical 
Research-Atmospheres, 114(17), doi: 10.1029/2009jd011702. 

Lammel M, Kroy K. 2017. Analytical mesoscale modeling of aeolian sand transport. Physical Review, E96: 052906, doi: 
10.1103/PhysRevE.96.052906. 

Li S, Li C, Yao D, et al. 2020. Wind tunnel experiments for dynamic modeling and analysis of motion trajectories of wind-blown 
sands. The European Physical Journal, E43: 22, doi: 10.1140/epje/i2020-11945-0. 

Ling Y, Wu Z. 1980. Experimention on the dynamic photography of the movement of sand-driving wind. Acta Geographica 
Sinica, 35: 174-181. (in Chinese) 

Lorenz R D, Zimbelman J R. 2014. Mechanics of sand transport. In: Lorenz R D, Zimbelman J R. Dune Worlds How Windblown 
Sand Shapes Planetary Landscapes. Berlin: Springer, 39-54. 

Pahtz T, Omeradžić A, Carneiro M V, et al. 2015. Discrete Element Method simulations of the saturation of aeolian sand transport. 
Physics, 42(6): 2063-2070. 

Reichardt H. 1941. About a new theory of free turbulence. Journal of Applied Mathematics & Mechanics, 21: 257—264. (in 
German) 

Shao Y, Raupach M R, Findlater P A. 1993. Effect of saltation bombardment on the entrainment of dust by wind. Journal of 
Geophysical Research: Atmospheres, 98(7): 12719-12726. 

Tan L, An Z, Zhang K, et al. 2020. Intermittent aeolian saltation over a gobi surface: threshold, saltation layer height, and high- 
frequency variability. Journal of Geophysical Research: Earth Surface, 125(1): e2019JF005329, doi: 10.1029/2019jf005329. 

Williams S H, Leer J A. 1995. Aeolian saltation transport rate: an example of the effect of sediment supply. Journal of Arid 
Environments, 30(2): 153-160. 

Yin Y. 1989. Study on sand drift in strong wind region in gravel desert. Journal of Desert Research, 9: 27-36. (in Chinese) 

Zhang K C, QuJ J, ZuR P, et al. 2008. Characteristics of wind-blown sand on Gobi/mobile sand surface. Environmental Geology, 
54: 411-416. 

Zheng X, He L, Wu J. 2004. Vertical profiles of mass flux for windblown sand movement at steady state. Journal of Geophysical 
Research: Solid Earth, 109(1), doi: 10.1029/2003jb002656. 

Zou X, Hao Q, Zhang C, et al. 1999. Parameter analysis on the saltation trajectory of wind blown sand. Chinese Science Bulletin, 
44: 1084-1088. (in Chinese) 

Zou X, Cheng H, Zhang C, et al. 2007. Effects of the Magnus and Saffman forces on the saltation trajectories of sand grain. 
Geomorphology, 90(1—2): 11-22. 


